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DISASTER++ is a C++ class library for the calculation of (1+1) and (2+l)-jet-like quantities in 
deeply inelastic lepton-nucleon scattering for one-photon exchange in next-to-leading-order 
QCD perturbation theory. The calculation is based on the subtraction formalism. The 
^v^j . user has access to an event record such that an arbitrary set of infrared-safe observables 

^ \ can be calculated in a single run. Compared to other existing universal programs, the 

full dependence on the number of flavours and on the rcnormalization and factorization 
^vq i scales is made explicit. An interface class providing a simple interface from C++ to existing 

FORTRAN programs is available. In a preliminary study DISASTER++ is compared to two other 
' programs for various bins of the lepton variables xb and y, where a particular emphasis is 

■ put on different behaviours for £ — > 1 of the parton densities /(C)- We find good agreement of 

DISASTER++ and DISENT (Version 0.1). The comparison of DISASTER++ and MEPJET (Version 
2.0) leads to several discrepancies. 
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1 Introduction 



For studies of the hadronic final state in high-energy collisions, versatile programs for the calculation of 
QCD corrections are required. The extraction of scale-dependent physical quantities such as the running 
strong coupling constant a s (Q 2 ) and parton densities fi (£,,Q 2 ) requires precise predictions in next-to- 
leading order of QCD perturbation theory. At the electron-proton collider HERA at DESY in Hamburg, 
the strong coupling constant has been measured via jet rates 0,[2|]. There is also a direct fit of the gluon 
density f g (£, Q 2 ) |§ based on a Mellin transform method |Q,|5|. Calculations for jet cross sections in deeply 
inelastic scattering for the case of the modified JADE scheme have been performed [f|-|To|lP| and implemented 
in the two programs PRO JET pi and D IS JET [Oj. 



In the meantime, calculations for arbitrary infrared-safe observables in deeply inlastic scattering have 



become available |13, 14| . In the last few years, the technology for the calculation of QCD corrections in 
next-to-leading order has developed considerably. The main problem in higher-order QCD calculations is 
the occurence of severe infrared singularities (they ultimately cancel for infrared-safe observables, or are ab- 
sorbed into process- independent, physical distribution functions such as parton densities and fragmentation 
functions). There are explicit algorithms available which permit the calculation to be done in a "universal" 
way: the infrared singularities are subtracted such that arbitrary infrared-safe observables can be calculated 
numerically. In principle, all existing algorithms are variations on a common theme, namely the interplay 
of the factorization theorems of perturbative QCD and the infrared-safety of the observables under con- 
sideration. There are two different ways to achieve the separation of divergent and finite contributions: 



the phase-space-slicing method and the subtraction method lq| . Both methods have their merits and 
drawbacks. 



(a) The phase-space-slicing method relies on a separation of singular phase-space regions from non-singular 
ones by means of a small slicing parameter s — > 0. The divergent parts are evaluated under the as- 
sumption that terms of C(s(log s) n ) can be dropped. The analytically evaluated phase-space integrals 
yield terms of the form (log s) m , which explicitly cancel against equivalent terms of opposite sign from 
a numerically performed phase-space integration. The simplicity of this scheme is obvious. The main 
problem is the residual dependence on the technical cut parameter s (in practice the limit s — > is not 
checked for every observable, but it is assumed that a fixed small value will be sufficient). Moreover, 
the numerical cancellation of the logarithmic terms by means of a Monte-Carlo integration is delicate. 
There is a calculational scheme available for the determination of the explicit phase space integrals 
over the singular regions [17]. For initial and final-state hadrons this scheme moreover requires the 



introduction of so-called crossing functions ]18[], to be evaluated for every parton density parametriza- 
tion. For deeply-inelastic lepton-nucleon scattering, an implementation of this calculational scheme is 



provided by Mirkes and Zeppenfeld in MEPJET [19]. 



(/3) The subtraction method is technically more involved, since the infrared singularities are cancelled point- 
by-point in phase space. The subtraction terms have, owing to the factorization theorems of pertur- 
bative QCD, a simple form. The problem is to arrange the subtractions in such a way that in the 
numerical evaluation no spurious singularities appear. A general framework, using a specific phase 
space mapping besides the factorization theorems, is given by Catani and Seymour in Ref. f20||, and 



implemented in D I SENT pi. 



i 



In these calculations, terms of the form clog c, c being the jet cut, have been neglected. This implies in particular a certain 
insensitivity with respect to the jet recombination scheme. The set-up of the calculations |], |?|,^ is such that a jet consisting 
of two partons is always mapped onto a massless jet. Therefore the jet definition scheme which is used on experimental data 
should be a "massless" scheme (this excludes, for example, the E-scheme). The variation of jet cross sections within the possible 
massless schemes cannot be modelled by that calculation. 
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The approach of the present paper is to use a generalized partial fractions formula to separate the 
singularities |2^| . The method is briefly explained in Section ^. We will describe in some detail the 
implementation DISASTER+-^ in the form of a C++ class library. 

There are two reasons for a new calculation, (a) The existing programs have the restriction that the 
number of flavours is fixed (Nf = 5 in the case of MEPJET and Nf fixed, but arbitrary for DISENT). 
For studies of the scale-dependence it is necessary to have a variable number of flavours, in order 
to be consistent with the scale evolution of the strong coupling constant and the parton densities. 
DISASTER++ makes the Nf dependence explicit in the "user routine" on an event-by-event basis, and 
thus results for several renormalization and factorization scales can be calculated simultaneously, 
(b) DISASTER++ is already set up such that the extension to one-particle-inclusive processes will be 
possible without the necessity of re-coding the contributions which are already present for the jet-type 
observables. This option will be made available in future versions of the program, as soon as the 
remaining contributions for one-particle-inclusive processes are implemented. 

The outline of this paper is as follows. In Section ^ we briefly review the algorithm employed in the 
present calculation. In Section ^ the FORTRAN interface to the C++ class library is described. Some remarks 
concerning the installation of the package are made in Section |j. A comparison of the available universal 
programs DISASTER++ (Version 1.0), MEPJET (Version 2.0) and DISENT (Version 0.1) is presented in Section ||. 
In a previous version of this paper, we have drawn the conclusion that we find an overall, but not completely 
satisfactory agreement of DISASTER++ and MEPJET, and that there are large deviations when comparing 
DISASTER++ and DISENT. One of the purposes of this paper is to present the results of a comparison of 
DISASTER++ and a new, corrected version (0.1) of DISENT. We now find good agreement of the two programs. 
We also give a few more results for MEPJET, in particular for the dependence on the technical cut s. It turns 
out that even for very small values of s we do not achieve agreement with the DISASTER++ / DISENT results 
for several cases under consideration^. The paper closes with a summary. The contents of this paper are 
mainly technical. The details of the calculation and phenomenological applications will be described in a 
forthcoming publication. 



2 The Algorithm 

The calculation is based on the subtraction method. A simple example to illustrate this method in general, 
and a comparison with the phase-space-slicing method, is given in Ref. pq| . For a more detailed exposition 
of the contents of this section, see Ref. 

The subtraction method is one of the solutions for the problem of how to calculate numerically infrared- 
safe observables without having to modify the calculation for every observable under consideration. In QCD 
calculations, infrared singularities cancel for sufficiently inclusive observables. The factorization theorems of 



perturbative QCD (see Ref. [26] and references therein) together with the infrared-safety of the observable 
under consideration guarantee that the structure of the limit of the convolution of the parton-level cross 
section with the observable in soft and collinear regions of phase space is well-defined and factorizes in the 



2 This is an acronym for "Deeply Inelastic Scattering: All Subtraction Through Evaluated Residues" . 

3 In a very recent paper E3L E. Mirkes quotes the results of the comparison of the three programs as performed in the 
previous version of this paper [M as resulting in a "so far satisfactory agreement" . This is a misquotation. The formulation in 
Ref. |24j] was that for MEPJET and DISASTER++ we find an "overall, though not completely satisfactory agreement", and that the 
results of DISENT (Version 0.0) "differ considerably" . Moreover, in the summary of Ref. p4| we mention that a few deviations of 
MEPJET and DISASTER++ are present. We wish to stress that there is a certain semantic gap between the expression "satisfactory 
agreement" and the results just quoted. 
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form of a product of a kernel and the Born term. This property allows, for the real corrections, the definition 
of a subtraction term for every phase-space point. Formally: 



soft/coll, limit , 



E 

A 



J dPS iA k A J dPS^- 1 ) 
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where a is the parton-level cross section, O is the infrared-safe observable, k A is a singular kernel, and t a is 
the non-singular part of the product a O. The index A runs over all possible soft, collinear and simultaneously 
soft and collinear singularities of a. The first integral is finite and can be calculated numerically. The second 
integral contains all infrared singularities. The term in the square bracket has a simple structure because 
of the factorization theorems of QCD, and the one-particle integral over the kernel k A and the factorization 
contribution from the term in the square brackets can be performed easily. This subtraction formula works 
only if the subtraction terms do not introduce spurious singularities for the individual terms that eventually 
cancel in the sum. This is achieved by a separation of all singularities by means of a general partial fractions 
formula 
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where the sum runs over all n\ permutations of n objects. 

In DISASTER++, the processes for (1+1) and (2+l)-jet production for one-photon exchange are imple- 
mented. The program itself, however, is set up in a much more general way. The implemented subtraction 
procedure can handle arbitrary number of final-state partons, and zero or one incoming partons (an ex- 
tension to two incoming partons is possible). The C++ class library is intended to provide a very general 
framework for next-to-leading-order QCD calculations for arbitrary infrared-safe observables. Of course, the 
explicit matrix elements (Born terms, virtual corrections and factorized real corrections) have to be provided 
for every additional process to be included. 



3 Program Structure 

We now describe the FORTRAN interface to the C++ class library. The C++ user interface will be documented 
in a forthcoming extension of this manual. 

To set the stage, let us first introduce some terminology. The user has to provide several subroutines 
which are called by DISASTER++ for every generated event. Each event e n , n = 1 . . . N consists of a set of 
phase spaces V nr , r = 1 . . . R n , and a set of contributions C n i, i = 1 . . . L n . Phase spaces V provide a set of 
four-vectors of initial and final-state particles, which are used to calculate observables 0{V). Contributions 
C n i consist of a list of weights w n ij, j = 1 . . . K n i (here: K n i = 11) which have to be multiplied by certain 
flavour factors F n ij. Every contribution C n i has an associated phase space V n r ni ', it is generally the case 
that particular phase spaces are used for different contributions. Flavour factors are products of parton 
densities, quark charges, powers of the strong coupling constant, and powers of the electromagnetic coupling 
constant. 

The expectation value (O) of a particular observable is given by the following sum: 

<°> = E E °(Pnr ni ) E WnijFmj. (3) 
n=l i=l j=l 

The first sum is the main loop of the Monte Carlo integration. 

The user has to provide a subroutine userl and a function user2. The subroutine userl (iaction) is called 
from DISASTER++ in the following cases: 
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iaction=l: 

after start-up of DISASTER++ 

iaction=2: 

before the end of DISASTER++ 

iaction=3: 

before the start of the grid-definition run of the adaptive Monte-Carlo routine, or before the final run 
of the adaptive integration, in case that there is no grid-definition run 

iaction=4: 

before the start of the final run of the adaptive Monte-Carlo routine 
iaction=5: 

after the final run of the adaptive Monte-Carlo routine 
iaction=6: 

once for every event (to initialize intermediate weight sums, etc.) 
iaction=7: 

signals that the event has to be dropped for technical reasons 

The function user2(. . .) is called from DISASTER++ after an event has been constructed. It has the 
following arguments (in an obvious notation): 

double precision function 



& user2( 

& integer nr, 

& integer nl, 

& double precision fvect(0. .3, -10.. 10, 1..30), 

& integer npartons (1 . . 30) , 

& double precision xb(1..30), 

& double precision q2(1..30), 

& double precision xi(1..30), 

& double precision weight (1 .. 11 , 1..50), 

& integer irps(1..50), 

& integer ialphas (1 . . 50) , 

& integer ialphaem(l . . 50) , 

& integer lognf(1..50) 

& ) 



Here nr stands for R n , nl stands for L n , f vect (mu, iparticle , ir) is the imr component of the four- 
vector of the particle with label iparticle (mu=0 corresponds to the energy component) in units of [GeV] 
in the Breit frame for the phase space ir; npartons (ir) is the number of final-state partons, q2(ir) is the 
value of Q 2 , and xi(ir) is the momentum fraction of the incident parton. The particle labels iparticle 
are given by 



iparticle=-8 : proton remnant 
iparticle=-7 : incident proton 
iparticle=-5 : outgoing electron 
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iparticle=-4 : incident electron 

iparticle=-l : incident parton 

iparticle=0. . (npartons-1) : outgoing partons 

The array weight (j , i) is a list of the weights for contribution i in units of [pb], irps(i) gives the 
index of the phase space for this particular contribution, ialphas (i) and ialphaem(i) are the powers of the 
strong and electromagnetic coupling constant, respectively, and lognf (i) is an index that specifies whether 
the weights have to be multiplied by a factor A consisting of a product of a logarithm of a scale and/or a 
factor of Nf. 

lognf =0: A = 1 
lognf =1: A = \n(n 2 r /Q 2 ) 
lognf =2: A = N f In (tf/Q 2 ) 
lognf =3: A = In (p 2 f /Q 2 ) 
lognf =4: A = N f In (pj/Q 2 ) 

Here p r and p,f are the renormalization and factorization scales, respectively. The total flavour factor for 
contribution i is given by 

F» 3 = A as ial P hasW o ial P ta « ft -, (4) 

where the quantity pij is a product of squares of quark charges Q a in units of e and parton densities. In 
particular: 

N f 

Pil= E Qlfa 

Q = l 

Nf 

Pi2 = E Qafa 
a=l 

Nf 2 
Pi3 — E Qafg 
a=l 

PiA = Pil 
Pi5 = Pi2 

Pi6 = Pil {Nf - 1) 

Pil = Pi2 {Nf - 1) 

Nf Nf 

Pi8 = E fa E Q% 
a=l P=l,/3^a 

N f Nf 

E fa E 

a=l p=l,p^a 



Pi9 = E fa E Q 2 p 



Nt N 



f 



PilO — E faQa E Qp 
a=l P=l,P^a 
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Nf N f 

Pill = Yl faQa J2 Q(3 
a=l P=l,P^ot 

The f a are parton densities evaluated for momentum fractions xi(irps(i)) and factorization scale pf, and 
fa stands for the parton density of the anti-flavour of the flavour a. The renormalization and factorization 
schemes are the MS scheme. The correction terms for the DIS factorization scheme will be implemented in 
the near future. 

We wish to note that the product of the weights, the flavour factors and the values of the observable is 
normalized in such a way that the sum yields the expectation value in units of [pb]. No additional factor 
such as 1/N, N being the total number of generated events, has to be applied in Eq. ^. 

Since phase spaces are used several times for different contributions, it is a good strategy to first evaluate 
the observable(s) under consideration for every phase space and to store the corresponding results. Then 
there should be the loop over the various contributions (the second sum). The innermost loop is the one 
over the flavour factors. 



The Monte Carlo integration itself employs the program VEGAS J27],|28|]. VEGAS is an adaptive multi- 
dimensional integration routine. Integrations proceed in two steps. The first step is an adaptation step in 
order to set up a grid in the integration variables which then steers the final integration step. The adaptation 
step itself refines the grid in a sequence of several iterations. VEGAS requires, as parameters, the number of 
Monte Carlo points to be used in the first and second step, respectively, and the number of iterations to 
refine the grid. In the framework of DISASTER++, VEGAS can be used in three different ways: 

• As an adaptive integration routine. The routine user2 should return a value. This value is handed 
over to VEGAS as the value of the integrand at the particular phase space point, and summed up. The 
final integral quoted by VEGAS is the sum of these weights for the final integration. This is the best 
choice if just one observable, for example a jet cross section, is to be evaluated. 

• As a routine that merely supplies random numbers for the events. If the number of iterations is set to 
zero, then VEGAS just performs the final integration run. The user is then responsible for the correct 
summation of the weights, and for the determination of the statistical error. It should be noted that, 
since all weights are available individually in the user routine, an arbitrary number of observables 
can be evaluated in a single run. In particular, since the dependence on the renormalization and 
factorization scales and on Nf is fully explicit, the study of the scale dependence of observables can 
be done in a very convenient way. For example, all plots from Ref. have been obtained in a single 
run of DISASTER++. 

• As a combination of the two preceeding alternatives. Here the adaptation steps are included. A 
"typical" infrared-safe observable, in the following called the adaptation variable, is evaluated, and its 
value is returned to VEGAS. This observable serves to optimize the distribution of points over phase 
space. A convenient observable of this kind is provided by DISASTER++ (see below). The "real" 
observables under consideration are evaluated as in the second alternative in the final integration step. 

DISASTER++ is initialized by a call of the subroutine disaster_ca() . It is recommended to end a DISASTER++ 
run by a call of the subroutine disaster_cb() in order to free dynamically allocated memory. 

Parameters can be set or commands be executed by means of three routines: 

disaster_ci(str, i) 

sets the integer parameter denoted by the character string str to the value i 
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disaster_cd(str , d) 

sets the double precision parameter denoted by the character string str to the value d 

disaster_cc (str) 

executes the command given by the character string str 

The following parameters are available (there are a few more to optimize the generation of the phase space 
points; they will be documented in forthcoming versions of this manual): 

ECM: 

the centre-of-mass energy in [GeV] 

LEPTON_INTEGRATION: 

1 : integration over x b and y 

XBMIN: 

minimum value of 

XBMAX : 

maximum value of xb 

YMIN: 

minimum value of y 

YMAX: 

maximum value of y 

QMIN: 

minimum value of Q 

QMAX: 

maximum value of Q 

WMIN: 

minimum value of W 

WMAX: 

maximum value of W 

PROCESS_INDEX : 

1 : leading order 

2 : next-to-leading order 

NUMBER_OF_FINAL_STATE_PARTONS_IN_BORN_TERM: 

1, 2, 3 for the process under consideration; 

1: (l+l)-jet-type observables (leading and next-to-leading order) 
2: (2+l)-jet-type observables (leading and next-to-leading order) 
3: (3+l)-jet-type observables (leading order only) 

POINTS: 

POINTS * (FACT_PREP + FACT_FINAL) is the number of generated points in the Monte Carlo integra- 
tion 

FACT_PREP : 

the number of points for the grid-definition run is given by FACT_PREP * POINTS 
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FACT_FINAL : 

the number of points for the final integration step is given by FACT_FINAL * POINTS 
RUN_MC : 

to start the Monte-Carlo integration 
A convenient adaptation observable can be evaluated by a call of the following function: 
double precision function disaster_cao( 



& integer ipdf _collection, 

& integer ipdf _parametrization, 

& integer ipdf_set, 

& integer i alphas _vari ant , 

& integer ialphas_order , 

& double precision dalphas_lambdaqcd4, 

& integer ialphaem_variant 

& ) 



The arguments of the function call are: 

ipdf collection: 

the collection of parton densities; 



1 : PDFLIB gg] 



ipdf _parametrizat ion: 

parametrization of parton densities (cf. PDFLIB) 

ipdf _set : 

set of parton densities (cf. PDFLIB) 

ialphas_variant : 

function which is used to evaluate the strong coupling constant; 

1 : running coupling a s (Q 2 ) with flavour thresholds at the single heavy quark masses 

ialphas_order : 

1 : one-loop formula 
2: two-loop formula 
for the running strong coupling constant 

dalphas_lambdaqcd4 : 

the QCD parameter Aq^ d for four flavours 

ialphaem_variant : 

function which is used to evaluate the electromagnetic coupling constant; 
1 : fine structure constant 
2: 1/137 

(an implementation of the running electromagnetic coupling constant is in preparation) 

To simplify the calculation of the flavour factors, a DISASTER++ routine can be called which returns the 
required coupling constants and the combinations of parton densities and quark charges: 
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subroutine disaster_cf f ( 



& 



integer ipdf _collection, 

integer ipdf _parametrization, 

integer ipdf_set, 

integer ialphas_variant , 

integer ialphas_order , 

double precision dalphas_lambdaqcd4, 

integer ialphaem_variant , 

integer nf , 

double precision ff actin(4) , 
double precision ffactout(13) 



& 



& 



& 



& 
& 



& 



& 



& 
& 
& 



) 



The arguments of the function call are the same as in the case of the routine disaster_cao (see above), 
except for the following: 

nf : 

the number of flavours Nf 

f fact in: 

input parameters; 

ffactin(l) : the momentum fraction variable £ 

f fact in (2) : the factorization scale in [GeV] (i.e. the scale argument of the parton densities) 

f fact in (3) : the renormalization scale in [GeV] (i.e. the scale argument of the running strong 
coupling constant) 

ff actin(4) : the scale argument of the running electromagnetic coupling constant 

f f actout : 

output parameters; 

ff actout (1 .. 11) : the quantities pn ■■■pm, 

ff actout (12) : the running strong coupling constant 

ff actout (13) : the electromagnetic coupling constant 

It is strongly recommended to use this routine, since it uses a cache that stores a few of the most recent 
values temporarily, such that the sums ptj and the parton densities do not have to be reevaluated. This 
routine is supplied for the convenience of the user. The weights and events generated by DISASTER++ do not 
depend on this routine. 

The description of the program structure just given may sound complicated. It is actually quite simple 
to use the program; an example for the calculation of the (2+l)-jet cross section for the JADE algorithm 
in the E-scheme is given in the files disaster_f . f and clust . f , as described in Section |j. 

4 Program Installation 

Source code: The source code of the class library is available on the World Wide Web: 

\protect\vrule widthOpt\protect\href {http : //wwwcn. cern. ch\string~graudenz/disaster .html}{hti 
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Files: The package consists of a number of files. To facilitate the installation, and to enable the C++ 
compiler to perform certain optimizations, the complete C++ part of the program is provided as one 
file onef ile_n. cc (the individual files are available on request). An example for the FORTRAN interface 
is given in the file disaster_f .f (calculation of the (2+1) jet cross section for the JADE algorithm 
in the E-scheme), together with a simple cluster routine in the file clust.f. The number of Monte 
Carlo events in the example is set to a tiny number (100) in order to terminate the program after a 
few seconds. Realistic values for the parameter POINTS are of the order of 10 6 . An example "make 
file" is given in makedisaster. 

Mixed Language Programming: DISASTER++ is mainly writen in the C++ programming language. The 
reason for the choice of this language are twofold: Object-oriented programming allows for programs 
that are easily maintained and extended^, and in high-energy physics there is a trend in the experimen- 
tal domain for a transition from FORTRAN to C++. Although the goal has been to write a self-contained 
C++ package, a few parts of the program are still coded in FORTRAN. Moreover, the standard parton 
density parametrizations are only available as FORTRAN libraries. This means that the DISASTER++ 
package cannot be run as a stand-alone C++ program. In most cases, users may wish to interface 
the program to their existing FORTRAN routines. An elegant and machine-independent way for mixed 
language programming for the case of C, C++ and FORTRAN is supported by the cfortran.h package 
described in Ref. [Q. For every FORTRAN routine to be called by a C++ method, an extern "C" routine 
has to be defined as an interface, and vice versa. The explicit calls are then generated by means of 
macros from cfortran.h. The most convenient way is, after compilation, to link the FORTRAN and 
C++ parts via the standard 



f77 -o disaster onefile_n.o ... 



command^, such that the FORTRAN part supplies the entry point. The required C++ libraries have to 
be stated explicitly via the -L and -1 options. The library paths can be obtained by compiling and 
linking a trivial program hw. cc of the type 

#include <stdio.h> 

mainO { printf ("Hello world! \n"); } 



with 



gcc -v hw. cc 



(for the GNU C++ compiler). An example for the required libraries can be found in the prototype "make 
file" makedisaster. Some machine-specific information is mentioned in the manual of cfortran.h. 

In the DISASTER++ package, the explicit FORTRAN interface, as described in Section ||, is already pro- 
vided. Thus from the outside the C++ kernel is transparent and hidden behind FORTRAN subroutines 
and functions. 

Template instantiation: In DISASTER++, heavy use is made of templates. At present, there is not yet 
a universally accepted scheme for template instantiations. The solution adopted here is the explicit 
instantiation of all templates. This requires that the compiler itself does not instantiate templates 
automatically. This is achieved for the GNU compiler by means of the switch 



-f no-external-templates 

4 It could even be said that object-oriented programming is a kind of applied ontology: the central categories of this approach 
are given by objects and methods that define their relationships. 
5 The procedure is described here for the UNIX operating system. 
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Output files: There is a small problem with the output from the C++ and FORTRAN parts of DISASTER++. 
It seems to be the case that generally C++ (FILE* stdout and ostream cout) and FORTRAN (UNIT=6) 
keep different file buffers. This is no problem when the output is written to a terminal, since then the 
file buffers are immediately flushed after each line- feed character. When writing to a file (as is usually 
the case for batch jobs), the file buffers are not immediately flushed, and this leads to the problem 
that the output on the file is mixed in non-chronological order. This problem will be solved by the 
introduction of a particular stream class which hands over the output to a FORTRAN routine. 

Miscellaneous: DISASTER++ employs the ANSI C signal facility to catch interrupts caused by floating 
point arithmetic. If the signal SIGFPE is raised, a flag in DISASTER++ is set, which eventually leads to 
the requirement that the event has to be dropped (via a call of userl(7)). Similarly a non-zero value 
of the variable errno of the ANSI C errno facility is treated. The signal handler is also active when 
the user routine is executed, which leads to the effect that in the case of a floating point exception 
the program does not crash, but continues under the assumption that the event has been dropped. 
Forthcoming version of DISASTER++ will make a flag available that can be used to access the status 
of the signal handler in the user routines. Moreover, it is checked whether the weight returned to 
DISASTER++ via user2 fulfills the criterion for IEEE NaN ("not a number"). If this is the case, it is also 
requested that the event be dropped. 

5 Comparison of Available Programs 

In this section, we compare the three available programs MEPJET (Version 2.0)0, DISENT (Version 0.1)0 and 
DISASTER++ (Version 1.0) numerically for various bins of xb and y as defined in Table [l], and for various 
choices of the parton density parametrizations. 





0.01 < y < 0.03 


0.03 < y < 0.1 


O.K y < 0.3 


0.005 < x B < 0.01 


Bin 1 (Q 2 > 4.6 GeV 2 ) 


Bin 2 (Q 2 > 13.5 GeV 2 ) 


Bin 3 (Q 2 > 45.0 GeV 2 ) 


0.05 < x B < 0.1 


Bin 4 (Q 2 > 45 GeV 2 ) 


Bin 5 (Q 2 > 135 GeV 2 ) 


Bin 6 (Q 2 > 450 GeV 2 ) 


0.2 < x B < 0.4 


Bin 7 (Q 2 > 180 GeV 2 ) 


Bin 8 (Q 2 > 540 GeV 2 ) 


Bin 9 (Q 2 > 1800 GeV 2 ) 



Table 1: Bins in xb o,nd y. The values in parentheses give the resulting lower bounds on Q 2 . 

The centre-of-mass energy is set to 300 GeV. To facilitate the comparison, the strong coupling constant 
is set to a fixed value of a s = 0.1, and the number of flavours is set to Nf = 5, even below the bottom 
threshold (Nf = 5 is hard-wired into MEPJET). The electromagnetic coupling constant is chosen to be 
a = 1/137 (the value is hard-wired into DISENT, but this could be changed trivially, in principle). The 
factorization- and renormalization schemes of the hard scattering cross sections are MS, and the factorization 
and renormalization scales \if and fi r , respectively, are set to Q. 

6 For the very-high statistics runs the default random number generator (generating a Sobol sequence of pseudo-random 
numbers) of MEPJET ran out of random numbers. We therefore had to modify the program such that it uses another generator 
which is also part of the MEPJET package. — The crossing functions for the "artificial" parton densities have been obtained by 
means of a modification of the program make_str_pdf 1 . f . 

' An earlier version of this paper |24j reported results of a comparison with DISENT Version 0.0. We found large discrepancies 
for some choices of the parton density parametrization. In the meantime, an error in DISENT has been fixed, and the results of 
DISENT and DISASTER++ are now in good agreement, see below. 
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The quantity under consideration is the (2+l)-jet cross section, shown in Tables 2-8 in Appendix |A|. 
For simplicity we consider the modified JADE clustering scheme with resolution criterion Sij <> cW 2 and 
the E recombination scheme, where Sij = (pi + pj) 2 , W is the total hadronic energy, and c = 0.02 is the jet 
resolution parameter. We require, in the laboratory frame (E e = 27.439 GeV, Ep = 820 GeV), a minimum 
transverse momentum of 1 GeV and a pseudo-rapidity of —3.5 < rj < 3.5 for all jets^j. 

The parton density parametrizations employed in the comparison are: 



(a) 


the MRSD' 


parton densities 


|31| 


(Table 2) 


(b) 


q(0 = (i - 


o 5 , 


9(0 = o 




(Table 3) 


(c) 


q(0 = o, 




9(0 = (1 - 


5 


(Table 4) 


(d) 


q(0 = (i - 


o 2 , 


9(0 = 




(Table 5) 


(e) 


q(0 = 0, 




9(0 = (1 - 


2 


(Table 6) 


(f) 


q(0 = (i - 


0. 


5(6 = o 




(Table 7) 


(g) 


q(0 = 0, 




9(0 = (1 - 





(Table 8) 



Here q(0 generically stands for valence and sea distributions^, and g(0 is the gluon distribution. We wish 
to point out that the comparison involving the "artificial" parton densities is not just of academic interest. 
On the contrary, for the extraction of, for instance, the gluon density from jet data it is convenient to replace 
the parton densities by simple functions with special properties (such as powers of the momentum fraction 
variable £ or functions of an orthonormal basis system), in order to achieve a fast fit. These functions 
usually do not have the shape of physical parton densities, in particular they do not have to fall off rapidly 
for £ — > 1. Moreover, next-to-leading-order calculations yield unique and well-defined results for the hard 
scattering cross sections to be convoluted with observables and parton densities. We employ the "artificial" 
parton densities also in order to have a stricter test of the hard scattering cross sections. 

The leading-order results of all three programs are in excellent agreement. The next-to- leading-order 
results of DISASTER++ and DISENT are in good agreement within about two to (sometimes) three standard 
deviations^ of the larger of the two errors quoted by the two programs. An exception is bin 7 for g(0 = 
(1 — £) 2 . A run °f DISENT with higher statistics yields a value of 0.1836 ± 0.0025, which is within two 
standard deviations of the DISASTER++ result, indicating that there was indeed a statistical fluctuation in 
the original result. 

The comparison of the next-to-leading-order results of MEP JET and DISASTER++ requires a more detailed 
discussion: 

8 These cuts in pr and r) are employed in order to facilitate event generation with MEP JET; the phase space generator 
implemented in that program is reminiscent of a generator for pp collider physics where pr and r\ cuts in the laboratory frame 
are a standard experimental procedure. It is thus complicated to generate events with MEP JET in the full phase space of the 
laboratory system, as usually required for eP scattering, where "natural" cuts in transverse momentum and pseudo-rapidity 
would be performed in the hadronic centre-of-mass frame or in the Breit frame. 

9 This means that «„(£), d v (Q, ^s(£), d s (£), s s (£), c s (£), b s (£) have been set to 

10 We wish to note that the error estimates quoted by the programs are usually not rigorous estimates because of the 
non-Gaussian distribution of the Monte-Carlo weights. Therefore, in principle, it is not possible to infer probabilities for the 
consistency of data samples produced by two programs based on these estimates. A more precise, but in general unfeasible 
way to obtain an estimate of the Monte Carlo error would be to run the programs a number of times with different random 
number seeds and to analyze the spread of the quoted results around their central value. Such a study has recently been done 
by M. Seymour for DISENT with the result that the individual error estimates are quite reliable E^]. 
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• For the MRSD'_ parton densities, the results for bins 3-9 are compatible within about two standard 
deviations of the statistical error of the Monte-Carlo integrations. The results for bins 1 and 2 differ 
considerably. Runs with a smaller value of the internal MEP JET cut-off variable s, which is set by default 
to s = 0.1 GeV 2 , yield the following results for bin 1: 580.6±6.7pb (s = 0.01 GeV 2 ), 564.8±10.5 pb (s = 
0.001 GeV 2 ) and 575.4 ± 13.0 pb (s = 0.0001 GeV 2 ). The statistical error is increased for decreased s 
because the integration volume of the (3+1) parton contributions is extended into the singular domain. 
Because of the increased statistical error, we also performed a high-statistics runs with ~ 4 • 10 9 (!) 
Monte Carlo events of MEP JET for this bin. For s = 0.001 GeV 2 we obtain 576.3 ± 6.7 pb and for 
s = 0.0001 GeV 2 the result is 583.2 ± 7.4 pb. These values from MEP JET are compatible with the 
DISASTER++ and DISENT results^- 

• For the parton density parametrization (b) (quarks only, with a steeply falling distribution g(£) for 
£ — ► 1), DISASTER++ and MEP JET are in good agreement. 

• The results for parametrization (c) (steeply falling gluon parametrization) are in good agreement, 
except for bin 1. 

• For parametrization (d), DISASTER++ and MEP JET are in agreement except for bins 1 and 4. Runs 
with a smaller value of the MEP JET cut-off variable s yield the following results for bin 1: 59.6 ± 1.8 pb 
(s = 0.01 GeV 2 ), 56.7±5.8pb (s = 0.001 GeV 2 ) and 54.9± 10.4 pb (s = 0.0001 GeV 2 ). A high-statistics 
run (~ 4 • 10 9 events) of MEP JET for bin 1 with s = 0.0001 GeV 2 gives the cross section 60.0 ± 1.9 pb. 
Contrary to the observation in case (a), for small s we do not get agreement of the MEP JET result with 
the DISASTER++ / DISENT result of about 48-49 pb. 

• The MEP JET results for parametrization (e) (g(S,) = (1 — £) 2 ) deviate considerably from the DISASTER++ 
results in bins 1, 2, 4 and 7. 

• For parametrization (f), DISASTER++ and MEP JET are incompatible for bins 1, 2, 4, 6 and 7. 

• For parametrization (g), MEP JET and DISASTER++ are compatible in bins 3, 5, 8 and 9 only. A high- 
statistics run (~ 4 • 10 9 events) of MEP JET for bin 4 with s = 0.0001 GeV 2 yields the cross section 
1.29 ± 0.02 pb. This value is different from the result for s = 0.1 GeV 2 , but still inconsistent with the 
DISASTER++ / DISENT result of about 0.69 pb. 

The overall picture is thus: Out of the three programs, DISASTER++ and DISENT (Version 0.1) are in good 
agreement within about two, sometimes three standard deviations of the quoted integration errors, both for 
"physical" and "artificial" parton densities. This agreement is very encouraging, but not yet perfect, and 
much more detailed studies involving different sets of observables and differential distributions are required. 
For the two programs, a direct comparison of the "jet structure functions" should also be feasible. 

For several bins, in particular for the "artificial" parton distribution functions, the MEP JET results for 
the default setting of the internal parameters deviate considerably from the DISASTER++ and DISENT results. 
For one particular bin studied in more detail for the MRSD'_ parton densities, the discrepancy disappears 
in the case of an extremely small internal technical cut s of MEP JET, for a substantial increase of the number 
of generated events to obtain a meaningful Monte Carlo error. A few MEP JET results employing "artificial" 
parton densities have been studied in more detail. We observed that in these cases a reduction of the s 
parameter does not lead to an improvement of the situation. For lack of computer time, we could not study 

11 These results underline that, for the phase space slicing method, results generally have to be validated ex post by a cross- 
check with a smaller technical cut s and much higher statistics. It may be argued that there are jet algorithms (the kr algorithm, 
for example) which show a better convergence for s — * 0. However, the point here is that one does not know in advance whether 
this is the case for the observable under consideration. — In Ref. [^3| we find the statement that s-independence in MEP JET is 
achieved for s = 0.1 GeV 2 . Our study shows that this is generally not the case, and that extremely small values of s, possibly 
of the order of s — 0.0001 GeV 2 , might be necessary. 
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all bins with a smaller s cut. The overall situation is thus still inconclusive and unclear. An independent 
cross check of the MEPJET results, in particular of those using the implementation of the crossing functions 
for the "artificial" parton densities, is highly desirable. 

6 Miscellaneous 

• If you intend to install and use DISASTER++, please send me a short e-mail message, and I will put 
your name on a mailing list so that I can inform you when there is a new version of the package. 

• Suggestions for improvements and bug reports are welcome. 

• In case that there are problems with the installation of the program, please send me an e-mail message. 

7 Summary 

We have presented the C++ class library DISASTER++ for the calculation of (1+1) and (2+l)-jet type ob- 
servables in deeply inelastic scattering. The program is based on the subtraction formalism and thus does 
not require a technical cut-off for the separation of the infrared-singular from the infrared-finite phase-space 
regions. A FORTRAN interface to the C++ class library is available. DISASTER++ is actually intended to be a 
general object-oriented framework for next-to-leading order QCD calculations. In particular, the subtraction 
formalism is implemented in a very general way. 

We have performed a comparison of the three available programs MEPJET, DISENT and DISASTER++ over 
a wide range of the parameters for the lepton phase space. We find good agreement of DISASTER++ and the 
Catani-Seymour program DISENT (Version 0.1). The comparison of DISASTER++ and the Mirkes-Zeppenfeld 
program MEPJET (for the MEPJET default parameters) leads to several discrepancies, both for physical and for 
"artificial" parton densities. For the MRSD'_ parton densities a reduction of the internal MEPJET phase-space 
slicing cut-off variable s, the number of Monte Carlo events kept fixed, leads to a certain improvement of 
the central values of the results, accompanied by a substantially increased statistical error and fluctuating 
central values. A considerable increase of the number of generated events (up to of the order of several 
billion events) eventually leads to an agreement of the MEPJET results with the DISASTER++ / DISENT results 
for a particular bin of the lepton variables which has been studied in detail. For "artificial" parton densities 
and a selected set of bins of the lepton variables, a reduction of the internal cut s does not resolve the 
discrepancies. Other bins are not considered for the lack of computer time for very-high statistics runs. It 
should be stressed that the present study is still limited in scope. An independent cross check of the MEPJET 
results for the "artificial" parton densities has to be done until a firm conclusion can be reached. Moreover, 
the study has to be repeated for a wider range of observables and much higher Monte Carlo statistics. The 
s dependence of the MEPJET results should also be studied in more detail. 

Compared to the other two programs MEPJET and DISENT, DISASTER++ makes the full Nj dependence 
and the dependence on the renormalization and factorization scales available in the user routine. This is 
required for consistent studies of effects such as the scale dependence when the bottom threshold is crossed. 
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A Numerical Results 



This appendix contains the numerical results which are discussed in Section ||. The entries in the tables are 
the (2+l)-jet cross sections in units of [pb]. 





leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


D I SENT 


DISASTER++ 


MEPJET 


DISENT 


1 


402 1± 1 13 

\J — • J 1 J. • -L *J 


399 9± 53 

KJ \J iJ • i_/ — 1 — KjiKJKJ 


399 6± 11 

KJ f. / * J . \J — 1 — -L • J. 


585 0± 2 6 

KJKJKJ • KJ — 1 — j—:,\J 


564 1± 19 

KJ KJ I • _L. 1 J. • i / 


578. 4± 7.1 


2 


207 6± 59 

1 • KJ — 1 — \J • KJ IS 


207 5± 34 

i-i\J t ' 'J — 1 — V 7 . KJ i 


207 4± 15 

j-J KJ 1 ■ ~ — 1 — v/ « i- K J 


364 8± 15 

KJ KJ X • KJ — 1 — ± • KJ 


347 3 ± 2 4 

f J 11 1 • • J 1 ■ i 


361 1± 3 5 

kj \j ± ■ j — i — kj • J 


3 


60.0± 0.16 


59.9± 0.14 


59.9± 0.15 


119.1± 1.71 


118.0± 1.05 


120. 1± 0.94 


4 


82.9± 0.16 


82.9± 0.10 


82.6± 0.21 


98.1± 1.11 


95.1± 0.61 


95.4± 0.87 


5 


42. 9± 0.08 


42.9± 0.06 


42.6± 0.28 


55.3± 0.46 


54.4± 0.49 


54.9± 0.40 


6 


11. 9db 0.02 


11.9± 0.02 


11.9± 0.08 


17.5± 0.06 


16. 8± 0.22 


17.3± 0.13 


7 


9.60± 0.03 


9.58± 0.01 


9.59± 0.04 


12. 1± 0.50 


12.7± 0.07 


12. 3± 0.15 


8 


6.24± 0.01 


6.23± 0.01 


6.24± 0.02 


8.61± 0.12 


8.55± 0.15 


8.52± 0.08 


9 


1.78± 0.003 


1.78± 0.003 


1.78± 0.06 


2.65± 0.03 


2.57± 0.06 


2.63± 0.02 



Table 2: Comparison for MRSD'_ parton densities. 
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leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


D I SENT 


DISASTER++ 


MEPJET 


DISENT 


1 


36.2± 0.09 


36.3± 0.05 


36.3± 0.12 


39.1± 0.33 


40.9± 0.89 


38.2± 0.53 


2 


17.8± 0.04 


17.8± 0.03 


17.7± 0.05 


23. 2± 0.37 


22. 7± 0.41 


22. 6± 0.22 


3 


5.21± 0.01 


5.21± 0.01 


5.21± 0.02 


8.24± 0.22 


7.86± 0.12 


8.14± 0.06 


4 


27.3± 0.06 


27.3± 0.03 


27.2± 0.09 


28.0± 0.52 


29.2± 0.18 


30.0± 0.21 


5 


14.8± 0.03 


14.8± 0.02 


14.7± 0.04 


17.4± 0.29 


16.9± 0.10 


17.0± 0.11 


6 


4.33± 0.008 


4.32± 0.006 


4.31± 0.01 


5.62± 0.10 


5.44± 0.05 


5.54± 0.03 


7 


6.38± 0.02 


6.37± 0.01 


6.38± 0.03 


8.49± 0.17 


8.59± 0.10 


8.37± 0.11 


8 


4.44± 0.01 


4.43± 0.007 


4.44± 0.02 


6.11± 0.08 


6.05± 0.07 


6.07± 0.06 


9 


1.36± 0.002 


1.36± 0.002 


1.36± 0.05 


2.02± 0.02 


2.00± 0.05 


2.01± 0.01 



Table 3: Comparison for q(£) = (1 — 





leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


DISENT 


DISASTER++ 


MEPJET 


DISENT 


1 


4.89± 0.017 


4.89± 0.007 


4.87± 0.01 


5.38± 0.07 


6.03± 0.06 


5.22± 0.13 


2 


2.66± 0.009 


2.66± 0.007 


2.65± 0.007 


3.67± 0.08 


3.66± 0.04 


3.58± 0.05 


3 


0.825± 0.003 


0.826± 0.002 


0.826± 0.002 


1.44± 0.07 


1.37± 0.03 


1.39± 0.02 


4 


1.60± 0.005 


1.60± 0.003 


1.60± 0.003 


1.20± 0.05 


1.30± 0.01 


1.12± 0.04 


5 


0.904± 0.003 


0.900± 0.001 


0.899± 0.002 


0.833± 0.027 


0.801± 0.008 


0.764± 0.019 


6 


0.279± 0.001 


0.278± 0.001 


0.278± 0.001 


0.314± 0.007 


0.287± 0.004 


0.299± 0.006 


7 


0.130± 0.001 


0.131± 0.001 


0.130± 0.001 


0.119± 0.005 


0.118± 0.002 


0.110± 0.006 


8 


0.0981± 0.001 


0.0980± 0.001 


0.0981± 0.001 


0.105± 0.002 


0.096± 0.001 


0.099± 0.004 


9 


0.0313^3.0001 


0.0310± 0.001 


0.0313± 0.001 


0.0396± 0.001 


0.034± 0.001 


0.0386± 0.001 



Table 4: Comparison for = (1 — ^) 5 
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leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


D I SENT 


DISASTER++ 


MEPJET 


DISENT 


1 


46.1± 0.11 


46.2± 0.07 


46.2± 0.14 


49.4± 0.67 


58.8± 0.65 


47.8± 1.2 


2 


23.8± 0.05 


23.8± 0.09 


23.8± 0.07 


30.6± 0.33 


31. 4± 0.71 


29.0± 0.54 


3 


7.28± 0.02 


7.28± 0.02 


7.29± 0.02 


11.2± 0.21 


11.0± 0.24 


11.4± 0.14 


4 


42.4± 0.09 


42.3± 0.06 


42.3± 0.12 


38.4± 0.30 


41.9± 0.26 


38.4± 0.31 


5 


23.9± 0.04 


23.9± 0.03 


23.8± 0.06 


24.8± 0.46 


24.2± 0.19 


23.9± 0.16 


6 


7.31± 0.01 


7.30± 0.01 


7.27± 0.02 


8.11± 0.19 


8.04± 0.41 


8.24± 0.05 


7 


20.3± 0.05 


20.3± 0.08 


20.3± 0.08 


23.3± 0.64 


25.1± 0.18 


22.4± 0.24 


8 


15.4± 0.03 


15.4± 0.02 


15.4± 0.01 


18.6± 0.36 


18.3± 0.47 


18.4± 0.15 


9 


4.87± 0.01 


4.86± 0.01 


4.87± 0.04 


6.47± 0.08 


6.38± 0.07 


6.41± 0.05 



Table 5: Comparison for q(£) = (1 — 





leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


DISENT 


DISASTER++ 


MEPJET 


DISENT 


1 


6.24± 0.02 


6.22± 0.01 


6.21± 0.02 


6.73± 0.13 


8.94± 0.12 


6.67± 0.24 


2 


3.59± 0.01 


3.58± 0.01 


3.57± 0.01 


4.77± 0.06 


5.24± 0.09 


4.43± 0.11 


3 


1.18± 0.004 


1.18± 0.004 


1.18± 0.003 


1.93± 0.04 


1.89± 0.04 


1.86± 0.03 


4 


2.65± 0.007 


2.65± 0.003 


2.65± 0.006 


1.13± 0.03 


1.66± 0.02 


0.94± 0.07 


5 


1.62± 0.004 


1.61± 0.002 


1.61± 0.003 


1.04± 0.04 


1.09± 0.02 


0.993± 0.03 


6 


0.535± 0.001 


0.534± 0.001 


0.533± 0.001 


0.433± 0.018 


0.412± 0.009 


0.430± 0.010 


7 


0.452± 0.002 


0.452± 0.001 


0.451± 0.001 


0.221± 0.026 


0.292± 0.010 


0.129± 0.02 


8 


0.398± 0.001 


0.398± 0.001 


0.397± 0.001 


0.298± 0.01 


0.271± 0.005 


0.237± 0.01 


9 


0.136± 0.001 


0.135± 0.001 


0.135± 0.001 


0.130± 0.003 


0.109± 0.002 


0.120± 0.004 



Table 6: Comparison for = (1 — 
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leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


D I SENT 


DISASTER++ 


MEPJET 


DISENT 


1 


50.6± 0.12 


50.7± 0.13 


50.7± 0.15 


58.6± 1.29 


72.9± 1.56 


54.7± 2.1 


2 


27.1± 0.05 


27.1± 0.16 


27.0± 0.07 


36.4± 0.57 


40. 0± 0.84 


34.9± 1.0 


3 


8.51± 0.02 


8.51± 0.02 


8.52± 0.02 


13.8± 0.35 


13.3± 0.43 


13.9± 0.2 


4 


49.8± 0.10 


49.7± 0.05 


49.6± 0.14 


41.2± 0.55 


47.2± 0.91 


41.9± 0.38 


5 


29.0± 0.05 


29.0± 0.03 


28.8± 0.07 


27.3± 0.52 


28. 2± 0.42 


26.4± 0.19 


6 


9.09± 0.01 


9.07± 0.01 


9.04± 0.02 


9.58± 0.06 


9.16± 0.15 


9.54± 0.06 


7 


30.6± 0.08 


30.5± 0.04 


30.5± 0.12 


32.0± 0.34 


36.3± 0.59 


32.4± 0.52 


8 


24.3± 0.04 


24.3± 0.03 


24.3± 0.07 


27.6± 0.56 


28.4± 0.35 


27.6± 0.21 


9 


7.88± 0.01 


7.86± 0.01 


7.87± 0.02 


9.63± 0.21 


9.50± 0.15 


9.47± 0.06 




Table 7: Comparison for q(£) = (1 — £) 




leading order 


next-to-leading order 


bin 


DISASTER++ 


MEPJET 


D I SENT 


DISASTER++ 


MEPJET 


DISENT 


1 


6.84± 0.02 


6.84± 0.01 


6.82± 0.02 


8.20± 0.25 


11. 6± 0.14 


8.26± 0.45 


2 


4.09± 0.01 


4.07± 0.01 


4.07± 0.01 


5.70± 0.11 


6.69± 0.16 


5.68± 0.17 


3 


1.39± 0.004 


1.39± 0.005 


1.39± 0.003 


2.41± 0.07 


2.33± 0.05 


2.34± 0.05 


4 


3.19± 0.01 


3.19± 0.01 


3.19± 0.01 


0.686± 0.09 


1.65± 0.03 


0.691± 0.10 


5 


2.06± 0.005 


2.06± 0.002 


2.05± 0.003 


l.OOrb 0.08 


1.14± 0.03 


0.866± 0.05 


6 


0.711± 0.001 


0.710± 0.001 


0.709± 0.001 


0.500± 0.006 


0.471± 0.01 


0.442± 0.017 


7 


0.712± 0.003 


0.711± 0.001 


0.710± 0.002 


0.157± 0.026 


0.373± 0.008 


0.082± 0.038 


8 


0.692± 0.002 


0.690± 0.001 


0.690± 0.001 


0.411± 0.020 


0.408± 0.022 


0.340± 0.023 


9 


0.245± 0.001 


0.245± 0.001 


0.245± 0.001 


0.194± 0.012 


0.172± 0.007 


0.161± 0.008 



Table 8: Comparison for g(£) = (1 — ^) 
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